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Abstract 



\ A generic ansatz is introduced which provides families of exact solutions 

to the equations of motion and constraints for Null strings in Bianchi type I 
cosmological models. This is achieved irrespective of the form of the metric. 
Within classes of dilaton cosmologies a backreaction mapping relation is es- 
tablished where the null string leads to more or less anisotropic members of 
the family. The equations of motion and constraints for the generic model are 
D ! casted in first order form and integrated both analytically and numerically. 



I. Introduction 

In recent years string theory has emerged as the most promising candidate for the 
consistent quantization of gravity and a unified description of all the fundamental 
interactions. A consistent quantum theory of gravity is necessary part of a unified 
theory of all interactions, because pure gravity (a model containing only gravitons) 
cannot be a physical and realistic theory. This is the strongest motivation for the 
study of strings in curved spacetimes [1-3]. 

In this context strings moving in curved spacetimes have gained considerable at- 
tention, since they could provide clues to a proper generalization of the theory [2-3]. 

Among other notable features, string theory has become a theory that gives inter- 
esting answers towards other fields, such as the physics of Black Holes, Cosmology, 
Galaxy formation etc [2-4]. 

In the cosmological context inflationary models have become natural outcomes of 
the theory. Yet a treatment of anisotropy which is a basic feature of Cosmological 
spacetimes has not be considered in the framework of string theory. 
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Anisotropic Cosmological spacetimes are candidates for a description of the early 
Universe, because it is known that inhomogeneous and anisotropic features may have 
prevailed in the primary stages of the evolution of the Universe. 

Since present day observations suggest that the Universe is highly homogeneous 
and isotropic, it would be interesting to examine whether string theory could account 
for a possible altering of anisotropic features that may have prevailed in the primordial 
stages of the Universe. 

Generic Bianchi I type cosmological solutions can occur in the context of more 
general dynamical theories [5]. As it is known dilaton fields appear naturally at the 
low energy limit of string theory, coupled with Einstein-Maxwell fields [6-8]. Also, 
dilaton fields appear as a result of dimensional reduction, of the Kaluza-Klein type 
Lagrangian [9,10]. The presence of a homogeneous primordial magnetic fields implies 
that the cosmological model is necessary anisotropic [5]. Much work has been done 
on magnetic Bianchi cosmologies in order to establish their main properties [11-14]. 
The introduction of the dilaton fiels renders the system in a completely integrable 
form [5]. It would be interesting to consider the general effects of the introduction of 
null strings in the above context, particularly in view of the possible interaction of 
null strings with the dilaton and/or EM fields and the fact that dilaton fields appear 
in the low energy-limit of string theory. The purpose of this paper is to study and 
examine these classes of models. 

The null string is a collection of massless point particles moving along the geodesies 
of curved spacetime [15]. Schild [16] was the first to consider null or tensionless 
strings (see also [17]). Null strings can be considered as the zeroth order term, of the 
perturbation expansion with respect to the string tension [18-20]. The evolution of 
null strings in different curved spacetime manifolds has been conidered extensively 
[21-23]. 

This paper is organized as follows: In section II we present the dynamical equations 
of bosonic and null string theory in general Bianchi type I models and we specify the 
classes of Dilaton-Bianchi cosmological spacetimes under consideration. 
The generic procedure for obtaining classes of null-string solutions is explained in 
section III and applied to the specified cosmological solutions. 

In section IV the energy-momentum tensor of string configurations is studied and the 
effects of string tension on the anisotropic features of these cosmological spacetimes 
is examined through the backreaction process. 

The thermodynamical properties of the general class of solutions, through the ansatz 
introduced, are studied in section V. 

In section VI, without further specification, the equations of motion and constraints 
for Bianchi type I models in their first order form, are integrated both analytically 
and numerically. 
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II. Strings in Bianchi type I Models 



The action for a bosonic string in a curved-spacetime background is given (for a 
D-dimensional spacetime) by [1] 

S = ~ J dTda^hh aP {T,a)d a X M dpX N G M N{X) (1) 

where (M,N = 0, 1, (D — 1)) are spacetime indices, (a, (3 = 0, 1) are worldsheet 
indices and T = (2na')~ l is the string tension. 

Variation of the action with respect to the "fields" which are the string coordinates 
X M (r, a), gives the following equations of motion and the constraints in the conformal 
gauge {h a p{T, a) = exp[(f>(r, a)\q a p) where <f>(r, a) is an arbitrary function and r] a p = 
diag(-l,+l) [2] 

X M - (X M )" + Y% [X A X B - X' A X ,B ] = (2) 
G AB (X) [X A X B + X' A X' B ] = (3) 
G AB (X)X A X' B = (4) 

where the dot stands for d T and the prime stands for d a . 

The equations of motion and constraints for the null-strings are given by 

X M + T%X A X B = (5) 
G AB (X)X A X B = (6) 
G AB (X)X A X IB = (7) 

The action (1) is invariant under general reparametrizations a a = <7 a (<7 /3 ) whereas 
the null string equations of motion and constraints only under T\ = f(r,cr), <J\ = 

9i?) [15]- 

We consider the four-dimensional Bianchi type I model as the background space- 
time, where the string coordinates are denoted by t = t(r,cr) , X = X(r,a) , 
Y = Y{r,a),Z = Z{r,a). 

The metric tensor is given by 

ds 2 = -dt 2 + a 2 (t)dx 2 + b 2 (t)dy 2 + c 2 (t)dz 2 (8) 

where a,b,c , are differentiate functions of the cosmic time. The Bianchi I space- 
time represents an expanding universe since the volume element of any spacelike 
hypersurface of simultaneity is constantly increasing [24,26] 



v^G = v/-(3)G = a(t)b(t)c(t) (9) 
The equations of motion for X^ = X M (r, a), (// = 0, 1, 2, 3) are given 
t-t" = ad[-(X) 2 + (X') 2 ]+bb[-(Y) 2 + (Y') 2 } + 

+ cc[-{Z) 2 + {Z') 2 ] (10) 

d T [a 2 (t)X}=d a [a 2 (t)X') (11) 

d T [b 2 (t)Y] = d a [b 2 (t)Y'] (12) 

d T [c 2 (t)Z] = d a [c 2 (t)Z'] (13) 
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The constraints become 



(i) 2 + (t) 2 ] + a 2 {t) [(X) 2 + (X') 2 ] + 
+b 2 (t) [(F) 2 + [Y'f] + c 2 (t) [(Z) 2 + (Z') 2 ] = 
-it + a 2 (t)XX' + b 2 (t)YY' + c 2 (t)ZZ' = 



(14) 
(15) 



We will be concerned with the case (t) > and normalize our solutions so that 
(r > 0). 

The open string boundary conditions demand 

(X M )'(r, a = 0) = (X M )'(r, a = vr) = (16) 

while closed strings 

X M (r, a = 0) = X M (r, <r = 2tt) (17) 

Introducing light-cone coordinates, x ± = (t ± ff ) i the equations of motion become 

d+d-t + aa(<9+X)(<9_X) + bb(d + Y)(d^Y) + 

+ cc(<9 + Z)(<9_Z) = (18) 

ad+d^X + d[<9 + X<9_t + <9+t<9_X] = (19) 

bd + d_Y + &[<9 + r<9_t + <9+t<9_y] = (20) 

cd + o_z + c[a+za_t + a + t9_z] = o (21) 

and the constraints read 

- \d + tf ± (cLt) 2 ] + a 2 (t) [(d + X) 2 ± (cLX) 2 ] + 

+b 2 (t) [(d + Y) 2 ± (cLF) 2 ] + c 2 (t) [{d + Zf ± (cLZ) 2 ] = (22) 

Recently, cosmological solutions of the Bianchi type I, possessing dilaton and elec- 
tromagnetic fields have been considered [5]. The solutions consist of two families of 
particular solutions for the scale factors and the dilaton field. The action is given by 
( note however several misprints in [5] ) 



S = J V=g(R + Xg^Wv + ^exp(-2uj)F^F a/3 )d 4 x 
and the equations of motion read 



(23) 
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(27) 


rp(MATT) 


= -A(w # w )V - -g^ v uj^ a ) - 






- exp{-2u){F m F au - l -g^F aP F af3 ) 


(28) 



with uj the dilaton field. The first family of particular solutions is given by 



a(r) = A e X0T 



c r 



6(r) = B e yoT 
C e z ° T , u(t) = Wvr + Wo 



(29) 



where the constants satisfy \Wi + 2(x y + y z + zoXo) = with A the dila- 
ton coupling constant and r defined with respect to the cosmic time t as dt = 
a(r)b(r)c(T)dr . The second family is given by 



a{r) = A r' p , b(r)=B r p 
c (r) = C t p , u(t) = W lnr 



(30) 



with p = A/(A — 2), W = 2/ (A — 2) . We want the kinetic term of the dilaton field 
to be positive, i.e. A < which requires < p < 1 . 

There exist also three families of general solutions. The first one which we shall 
deal with, is given in terms of the initial expansion rates along the three directions 
(001,002,003) . Here Qo = 001002 + 002003 + ooioo 3 , r 2 = A 2 {1 + A 2 ) , A = — 2A 2 and 

S = A 2 LO\ . 



a{r) 
b(r) 
c(r) 



s A 2 
-— t + oo -In 



smhl—r 



A^e-^l^ 



C e {sA2/r2)T 
+ (s/r) 2 



[a^ t , 

• t / fir 

[a^ t , 

^ -r 



-1 A 4 /r 2 



-1 A A /r 



(31) 

(32) 

(33) 

(34) 
(35) 



III. The Generic Ansatz and Exact solutions 

We denote the target space coordinates of the string by t = t(r, a), X = X(t, a),Y = 
Y(r,a),Z = Z(r, a) , and introduce the ansatz, t = t[(j)],X = X [(/>], Y = Y[<j)],Z = 
Z[(f>] , for the string coordinates, with = 0(r, a) a function to be determined. Hence 
we postulate an explicit functional-form dependence for the string coordinates. 
The metric tensor functions of Eq (8) assume the form 

a 2 (t) = a 2 [t( ( f ) )} = a 2 (<j)) (36) 
b 2 (t) = b 2 [t(<j ) )] = b 2 ( ( j ) ) (37) 
c 2 (t) = c 2 [t(0)]=5 2 (0) (38) 
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Now we can define the functions £ = £[0(r, a)], r\ = r][<f>(T,cr)], ( = £[0(r, er)] , 
through the differential relations 



d£ = a((f))dX = a(</>)X t4> d<i> (39) 
dr] = b(<f>)dY = b((f))Y^d<f) (40) 
d( = ~c(<f))dZ =~c(<P)Z^d<p (41) 

The comma's with respect to <fi denote total differentiation of the indicated functions, 
these being considered as single- variable functions of it. The constraints (22), become 

& + Vl + & - tl = (42) 



Eqs (19)-(21) become 



+ 



a 



vAd+d-^) + V,ct><t>(d + <f))(d-<f)) + 
~b t 

+ [t 

+ 



t^O(9 + 0)(9_0) = o 

) + 0)(d_0) + 

V^(d + 0)(cL0) = O 

) + 0)(9_0) + 



Eq (18) yields 



tAd + d-<t>)+t M {d + <t>){d-<t>) + 



a 



+ 



b 



M 2 (9 + 0)(9_0) + 



c 



(O) 2 (9 + 0)(9_0) + 

(O) 2 (<9 + 0)(<9_0) = O 



We consider here three classes of solutions, corresponding to the choices 



I. 
II. 

III. 



(d + d^) = 

(d + cL0) = /<d + 0)(d_0) 



(43) 
(44) 
(45) 



(46) 



(47) 
(48) 

(49) 



for the function = 0(x ± ) • Here ^ are non-zero constants. The general 

solution of the first equation (wave equation) is given by 



(50) 
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while for Eq(48), 



= 0(X ± ) 



-1 



In 



H{x + )+G{x') 



(51) 



where H, G are arbitrary integration functions of the indicated arguments. 

I. The equations of motion (33)- (36) and the constraint (32) become a set of four 
second-order ordinary differential equations and a first order differential constraint, 
namely 



t.d>6 + 



a 



(C 



+ 



b 



M 2 + (O) 



o 



(52) 



and 



Eqs (53)- (55) integrate to 



V,<M> 

Cm 



01, t 



L a 
b 

c 



; 

a 



V,4> 



t A,<t> = 
t,4,ri,<t> = 
t ,<t>C,4> = 

t 2 -t 2 - 

c y /■ _ 
: "~T> s,<* — 



2i 

c 



(53) 
(54) 
(55) 

(56) 
(57) 



where c x ,c y ,c z are arbitrary constants, while substituting Eqs (57) into (52), and 
integrating once, results in the same equation to be satisfied with the case when one 
substitutes Eqs (57) into (56), which is 



[t,4 



a 2 it) b 2 (t) c 2 (t) 



(58) 



Therefore for every given set of scale factors a(t), b(t), c(t) one solves the ordinary 
differential equation (58) for t = t((f>) . Then through Eqs (57) one determines 
£ = £(0)) r 7 = ^(0)? C = C(0) ) with = 0(x ± ) satisfying Eq (47). Finally from Eqs 
(39)-(41) one obtains the explicit expressions for the string coordinates. 
II. Following the same line of reasoning, we obtain 



a b 
where c x ,c y , c z are arbitrary constants and 



K<f> 



(59) 



+ 



+ 



a 2 [t) b 2 (t) c 2 (t) 



-2n<p 



(60) 
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Now <fi = 0(x ± ) satisfies eq (II) or equivalently Eq (51). 
III. In the same manner, we obtain 



P — n = 2]L t = — 

™ Va' V '+ W ™ t»c 



where 



c x ,c y , c z are arbitrary constants and 



+ 



a 2 (t) b 2 (t) c 2 (t) 



(61) 



(62) 



Now <fi = 0(x ± ) satisfies Eq (49) after t = t(<f>) as the solution of Eq (62) is 
substituted. The interesting fact is that given any set of Bianchi-type Scale-factors, 
irrespective of the general dynamical theory that they come one proceeds directly to a 
string solution depending on two arbitrary integration functions and seven constants. 
The invariant string size is defined as the induced metric on the worldsheet [3,15] 



ds 2 = G AB {X)dX M dX N = G AB {X)X^X N A-dr 2 + da 2 ) 



M V N i 



(63) 



which vanishes due to the constraints (22). Therefore the worldsheet hypersurface 
is null and we have classes of string solutions that are null-string solutions. It is 
straightforward to verify that Eqs (5)-(7) are satisfied for the choices of the = 
, if we insert the small dimensionless parameter of the null-string perturbation 



expansion c 2 = 2AT (0 < c < 1) [19] in the terms that contain the (a) 
derivatives in Eqs (57)-(59) and let (c — > 0) . For example Eq (59) will become 

(0-cV) = M^A)((0) 2 -c 2 (0') 2 )- 

However, if we let the parameter (c) to be small (c< 1) we can always reconcile 
the fact that the solutions are null-string solutions (i.e. have null worldsheet manifold) 
with the fact that they also satisfy the full bosonic string equations of motion (2)- 
(4). Indeed, for this case while the solutions satisfy the full string equations (2)- 
(4), when expanded with respect to the parameter (c) they also satisfy in a very 
good approximation the null equations (5)-(7). This will permit us to use the string 
energy- momentum tensor derived from the action (1) and study the effects of the 
string tension on the anisotropic features of spacetime, at least in the first order case 
with respect to the perturbation parameter (c) . Thus these classes of solutions are 
more general than solutions that arise from Eqs (5)- (7). 

For the second family of particular solutions, Eqs (30), in the case of p — 
+ (1/2) and p = (1/6) the equations are integrable in closed form. One obtains 
the following, performing the integration for (I) (and in a completely similar way for 
the cases II, III, Eqs (48)-(49) ), 
1. For p=+(l/2) 



a(t) = A 
b(t) = B 



St 



2A B C 
3t 



-1/3 



= a r 1/3 



2A B Co 



1/3 



= b t 1/3 



(64) 
(65) 
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and for the solution 



2. For p=+(l/6) 



and for the solution 



c(t) 



Co 
2W 



3f 

2A B Co 
In ' 



1/3 



= c t 1/3 



2A B C 



^£,4/3 , 2 
2 l ~ ^2/2 



1/2 

^ = 

Z ,<t> = 

c 2 = 

^yz — 

t = 



+ 



t 4 / 7 + c; 



[cr/a 2 (0] 
[c,/6 2 (t)] 

(4/bl) + (cl/cl) 

t{4>) 



a(t) 


= A 


b(t) 


= B 


c(t) 


= C 




6W 
7 



7t 



SA B C 
7t 



6A0B0C0 
7t 



1-1/7 
1/7 



QA0B0C0 
■In ' 



GA0B0C0 



\ 1/2 (c 2 

S/2 



«0 



^+4/7 _ o r 2 
2 ^ * u yz 



"1(0 + 00) 



(66) 
(67) 

(68) 

(69) 
(70) 
(71) 
(72) 
(73) 

(74) 
(75) 
(76) 
(77) 

(78) 



X i4 , = [c x /a 2 (t)] 

Z A, 



t 



[c y /b 2 {t)\ 

[Cz/C 2 (t)} 
t(0) 



(79) 
(80) 
(81) 
(82) 



As an explicit demonstration of the integration procedure for the case p — (1/2) we 
have the complete solution for the choice <po = c yz , as follows. We insert Eq (68) 
into Eqs (69)- (71) implicitly, through Eqs (64)- (66) and obtain 



1/2 



4c, 



0+1 



3«q 



+ 
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(83) 
(84) 
(85) 



where satisfies the wave equation (47). 

However we will be concerned with the most general case Eq (39) and keep the 
parameter p unspecified. We have 



a(t) 
b(t) 
c(t) 




Co 
W 



i -p/(i+p) 



p/(i+p) 



-i p/(i+p) 



1+p 



(l+p)* 
AqB C 

" (l+p)t 



■In 



(86) 

(87) 

(88) 
(89) 



and 



f — t-(p+n+pft)/0-+p) 



t 4 P /(l +P ) +c 2 



1/2 



t 



[c y /m)\ 

[c z /t»c{t)\ 



(90) 

(91) 
(92) 
(93) 
(94) 



The cases that are integrable in closed form corresponding to the above two cases are 
precisely p — (— 1 + 2p)/(l +p), and p = (— 1 + 6p)/(l +p) , but we will retain the 
general form throughout. 

For the case of the solution (30), from Eq (28) we get (following the standard 
procedure of [5]) 



87rT f 



(MATT) 



A 



w 2 



F 2 



00 



2(l+p) 2 t 2 2B 2 C 2 Q 



8vrTi MATT ) = 



A 



w 2 



F 2 



2{l+p) 2 t 2 2B 2 C 2 



(l+p)t 

AqBqCq 

OTpji 



(W +2p) 



(W +2p)' 



(95) 



a 2 {t) (96) 



10 



a T (MATT) = 

yy 



8nT (MATT) 



A W 



F 2 



A Wl 



n 



2(l+p) 2 t 2 2BICI 




TiT5T(Wb+2p)- 



b\t) 



c\t) 



(97) 
(98) 



IV. The Energy-Momentum Tensor 



The energy-momentum tensor for the string is obtained from the functional derivative 
of the action (1) at the spacetime point X. We obtain [4] 

V=GT AB {X) = ~f drda(X A X B - X' A X' B )5 {A \X - X(t, a)) (99) 



where 

5^(X - X(t, a)) = [] KX M ~ X m (t, a)) (100) 

M=0 

It is useful to integrate the energy-momentum tensor in a three-dimensional volume 
that completely encompasses the string [3] 

I MN (X) = J ^/^GT MN {X)d 3 X (101) 

For the ansatz introduced, t = t[0(r, a)] , X — X[<j>], Y — Y[<j>], Z = Z[4>] we obtain 



rMN 



VO--2J 



T r dadt{r, a) 



X^X-UY-^)% =m 5[t-t^a)] 



(102) 



Finally using Eq (9), an effective energy- momentum tensor for the string is derived 



iMN 



(STRING) 



(X) 



Ana' a(t)b(t)c(t) 
[(0) 2 -(0') 2 



t_rh 



J dadcj) 







t=t[<t>] 
T,0-)} 



(103) 



We demonstrate the evaluation of this integral for the case that = 0(x ± ) satisfies 
the wave equation (47), for a closed string. We have 



J dadcf) 



'\21 



r2n 

= Jo ' 



da 



m 2 - (0 ) 

(0) 

Uf - (0?: 

(0) 



-5[<f> - 0(r, a)} = 



(104) 



t=t{4>,<t) 
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where the dependance of (r) is now in the variable (0) defined by the cosmic time 
t = t((f>) . For closed strings satisfying the wave equation we have 



= + 2qa t + 



n 



(105) 



with the reality condition 0* = 0_ n {n G Z) and a similar relation satisfied by the 
tilded operators. We can invert approximately this relation to obtain r ~ (<p/2qa ) 
and finally 



I = Anqa H Re[(j) n <j) n exp(— 2mr)] 

9 n>0 



(106) 



In the general case, one must evaluate 



= y dcrd0 



[( 



(0) 



-5[(f) - 0(r, <t)] = y d<rd*(r, 



[(«^) 2 - (0?] 



(0) 



lT ;_- _ -S[t-t(T,(T)] 

(107) 

where = 0(x ± ) satisfies Eq (39). This task is difficult to be achieved in its full 
generality. However we can estimate the dependence on the cosmic-time t . From 
Eq (49), neglecting the primed factors we obtain approximately [ln{4>) 2 }^ = [Int^]^ 
and since the fraction in Eq (107) is roughly proportional to (0) we have / oc . 
Therefore we obtain 



T(STRING)(t) (X 



4W a(t)b(t)c(t) 



M vNl 



t. 



(108) 



4>=<P(t) 



With the aid of Eq (108) we write the functional dependence of the energy- momentum 
tensor components of the string and the "matter" content of the model that corre- 
sponds to Eqs (86)-(94). 



T 00 nr 
1 (STRING) <x - 

1 (STRING) °^ 



t (p+^+pm) /(i+p)fp/0-+p) 

t p/(l+p) 



tl t/2 



rpOy 

1 (STRING) 



T (STRING) « t 3p /(1+p)f/i/2 

t (p+^+p^)/(l+p) t Zp/(l+p) 



1 



1 (STRING) 



rpXX 

1 (STRING) 0C 



rpzz 

1 (STRING) 



t 3„/2 /( c 2/ a 2 )t 4 P /(l +P ) + c ^ 



(109) 

(110) 
(111) 

(112) 
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= T, 



yz 

(STRING) 



oc 



rpxz 

1 (STRING) 



^(STRING) 0C 



t 5 P /(l+ P ) t 3 M /2 ^(4/^2)^/(1+^) + c : 

^(p+^+pm)/(i+p) 



.2 



t p/(l+ P )t3 M /2 /( c 2/ a 2) t 4p/(l+p) + C 2 



(113) 



(114) 



Correspondingly, for the "matter" content of the model, they are given by Eqs (95)- 
(98) with T xx = a~\t)T xx , T yy = h~\t)T yy , T zz = c~\t)T zz , using the fact that 



< p < 1 with A 



07ri (MATT) 



Q rpxx 

° nl (MATT) 



(2p/(p-l)) <0 



F 2 



(l-p)(l+ p )2 t 2 

1 pW 2 



A 2 (l-p)(l+p)^ 2 
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A) -So Co 



;i+p)t 



- (W +2p) 



(i+p) 



(115) 



+ 



anl (MATT) 
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AqBqCq 

JT+pji 



TTT5)(^'+2 P )- 



AqBqCq 



B 2 (1 - p)(l + p)H 2 



& n T*MATT) 



o r2^2 
1 



(l+p)f 

v wl 



TTT5)( H/ o+2p)- 



Aj-BqCq 



TT+pT 



(116) 



C 2 Q (1 - p)(l + p)H 2 



2B 2 Cl 



AqBqCq 



:i+p)t 



A)-SoC 



2p 



(117) 



2p 



(118) 



-oo 



For small values of the cosmic time (t — > 0) or for large values, (£ 
demanding functional identification, between the string energy-momentum density 
(00) component, Eq (109), and the EM part of the "matter" content, Eq (115), 
results in 

4(^o +p) 



// 



(l+p) 
4(W^o + 2p) 



(f-0) 
(t -> +oo) 



(119) 
(120) 



The limit of t — ^ is taken by neglecting the first term in the square root of Eq (109), 
while the limit of t -> +oo by neglecting the second term. The same argument, 
regarding the pair of dilaton and EM parts of the "matter" content in Eq (115), yields 
W — (1 — p) , therefore, 

4 



(l+p) 
= 4 (t 



(f-0) 
-> +oo) 



(121) 
(122) 
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For the second family of particular solutions Eqs (30), when the Universe is not highly 
anisotropic, we have (p ~ 0) and so the components of the string energy-momentum 
tensor Eqs (109)-(114) do not exceed the (00) component Eq (109), as t approaches 
zero. We can now form a sort of backreaction process, regarding the EM field F as 
a small perturbing factor that can be ommited. The coefficients of the dual Universes 
are related by 

p{1 - p) ±( T )c - W-rt (123) 

In Eq (123) T = (l/2na') is the string tension. Using Eq (109), we assume that 
the presence of the null-string acts as an additional new source, so that the primed 
term corresponds to another member of the family, Eq (115). Eq (95) has also been 
used and Wq = (1—p) is substituted. This relation is valid in the immediate vicinity 
of (0 ~ p < 1/3) . We obtain now the range of the cosmic-time where this coupling 
occurs. The fraction at the r.h.s. of Eq (123) is dimensionless and has a maximum 
equal to (1/8), at p = (1/3) . So the last term on the l.h.s. of Eq (123) has 
to be dimensionless and of the same order. Using Eq (90) we can therefore write 

Tt (4+p)/(l+p)^ _ A 

It is well known that the string tension can be taken to be of the order of the Planck 
energy scale [3] T Pl = (1/2W) ~ 10 38 (GW) 2 where the Planck mass is given 
by M pianck = (hc/G) 1 / 2 ~ 10 19 '(GeV/c 2 ) and correspondingly the characteristic 
time scale of quantum gravity is tgc — 10 _43 sec while for the GUT's we have 
tour — 10~ 30 sec . The gravitational constant is G ~ 10~ 30 (GeV)~ 2 . Since we have 
null strings we use in the action (1) the tension T to be much smaller than the one 
for the tensionful strings thus (T/Tpi) ~ 10~ K with («) a large integer. 

The quantity has the form of inverse angular frequency, and for closed strings 
it is proportional to the modes of oscillation that are excited. We can therefore set 
equal to a integer n, (n=l,2,...). So we obtain 



T 



TpiJ \tQG, 



(4+p)/(l+p) 

| ~ O.ln (124) 



Taking logarithms we obtain 

/4 + p" 



log w - — = k - 1 + log 10 n (125) 



\l+pj \t QG/ 

This formula provides the range of cosmic-time, with respect to the used value of 
the string tension, where the coupling occurs. As the value of the string tension is 
increased (k — > 0) or we have excitation of higher modes of the string (increasing n), 
the range approaches the Planckian regime. We conclude that when the string action 
enters with a negative sign, Eq (1) we get a transition to less anisotropic member 
of the family but also because of the presence of the dilaton field in this family of 
solutions it exhibits also a sort of chaotic behavior, in which case one can get to a 
more anisotropic member of the family (see Fig 1.). 
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V. The General Case 



We consider now the first family of general solutions Eq (31)-(35). These correspond 
to a Universe which is singular at t — > , where the energy density of th magnetic field 
and the scalar fields diverges and the spatial volume tends to zero [5]. Using Eq (12) 
of [5], it is not possible to obtain in closed form the dependance of the scale-factors 
on the cosmic-time t. (note that it should be written as dt = a{T)b{T)c{T)dr) . 
However we can examine the cases of the strong and weak coupling limits of the 
Dilaton field. These correspond to the limits A 2 — > +oo and A 2 — > respectively. 
We obtain 




1 

sinh(f3r) 
1 



(126) 
(127) 



For the strong coupling limit we get 



t = 



P 2 



;[Pcosh(Pr) — usinh{l3r)} + constant 



v = (2uj 1 + u 2 + uj 3 ) 



(128) 
(129) 



where the integration constant is chosen so that the cosmic and worldsheet times to 
assume zero value simultaneously. Now we can study the asymtotic behaviour of the 
solutions in two particular cases: 



r-0, t~-^2-^[e"-l] (130) 
Here we get (3 = (u^a^ + cj 2 ^3 + ^3^1) 2 + oj\ 



-In 



1 + 



(P 2 



A 



a ^ = sinh(pr) 
b{t) = B e^ 1+ ^ T sinh{(3r) 
c (t) = C e^ 1+UJ ' i)T sinh(pr) 



(131) 



The energy-momentum tensor for the string is given by Eq (108) together with Eqs 
(62) and (91)-(93). We now give a physical interpetation of the parameter (//) , in 
Eq (49), which is related to the thermodynamical interaction of the string with the 
geometry. From Eq (130), close to r ~ we have that r ~ [(f3 2 — v 2 )/isf3]t and the 
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scale factors become 

a(t) = 



A u 



(P 2 - v 2 )t 

(P 2 - ^ 



b(t) = B ^- — >-t 

?2 ,,2 



(0 2 - v 2 \ 

c(t) = C [P ' t (132) 

From Eq (108), inside the three dimensional volume that completely encompasses the 
string, V = a(t)b(t)c(t) = A B C (/3 2 - v 2 )tjv , using Eqs (101)-(103) we obtain the 
total energy E^ STR \t) of the string and the pressures as follows 



E (STR)^ 



r 2 c 2 r 2 

^x _|_ _|_ u z 

a 2 (t) 6 2 (t) c 2 (t) 



1/2 



Px 
Py 



1 



,-,2 



4W tn/ 2 Va 2 (t) 

1 c 2 
1 c y 

'Ana' t»/ 2 Vb 2 {t) 



1 c 2 

Pz = "W ^7W(0 (133) 

It is well known [25] (see also [26], p. 449) that if F^ STR \t) is the /ree energy of the 
string and /3 the string temperature, we have 

E (STR) = F (STR) + (l/fyS {STR> 

p k = -(dF^/dV)? (134) 

where 5'( s " rR ) = f3 2 [dF^ STR) / d(3) is the string entropy, in the three-dimensional 
volume. So if we consider isothermal motion of the string, then its temperature does 
not depend on time. Using Eq (133), and the fact that the volume (V) is proportional 
to the cosmic time t, we have 

oc 135 

dt 2 V K J 

so it is a measure of the entropy exchange between the string and the spacetime 
geometry. 

If we consider adiabatic motions of the three-dimensional volume that contains 
the string i.e. one allows for time dependance of the temperature /3 = /3(t) in order 
that there exists no entropy exchange then one concludes that ji is a measure for the 
adjustment of the string temperature that is required for no entropy exchange. Using 
again Eq (134) we demand that there exists no entropy exchange dS^ STR ^/dt = 
resuling in 

~w K (136) 
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VI. Dynamical Isotropization and Flatness by Null 
Strings 

In this Section, we integrate numerically the Einstein field equations which determine 
the evolution of a Bianchi Type I model, with matter content consisting solely of null- 
strings. A null-string may be considered as representing a collection of points moving 
independently along null geodesies [16]. The corresponding equations of motion and 
constraints are now given by Eqs. (5) - (7). In this case, Eq. (7) ensures that each of 
these points propagates in a direction perpendicular to the string. 

In general, it is expected that the introduction of a null-string structure in a 
cosmological model will react back on the curved background, modifying its dynamical 
characteristics. In what follows, we are interested in determining the way that this 
backreaction procedure affects on the evolution of the Universe. Accordingly, we need 
to calculate the components of the energy-momentum tensor attributed to a spacetime 
region due to the propagation of a null-string, in an otherwise vacuum cosmological 
model, with metric in the form of Eq. (8). These quantities are subsequently inserted 
into the r.h.s. of the corresponding field equations. The resulting dynamical system 
is highly non-linear and therefore, solutions can be obtained only through certain 
numerical techniques, where the concept of attractor plays an important role: If some 
special spacetime represents an attractor for a wide range of initial conditions, such 
a spacetime is naturally realized asymptotically. 

In a Bianchi Type I spacetime, the equations of motion for a null-string in the 
zeroth order approximation to c 2 = 2AT, are 

dec o _ dib o d/C o , \ 

t + a — x 2 + b — y 2 + c — z 2 = (137) 
at at at 

1 da ■ 
x + 2 - — tx = 
a dt 

1 db .. n 

» + 2 Zdt tv = 

1 dc ■ , 
z + 2 - — tz = (138) 
c dt 

where a dot denotes derivative with respect to the world-sheet coordinate r. The 
corresponding constraints [Eqs. (6) and (7)] are written in the form 

i 2 = a 2 x 2 + b 2 y 2 + c 2 z 2 

it' = a 2 xx +b 2 yy' + c 2 zz' (139) 

where the prime denotes differentiation with respect to a. Eqs. (138) are evaluated 
as follows 

x (a) 

x = — ~ 

a 2 
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b 2 



c 



.2 



(140) 



In this case, provided that t ^ 0, from Eq. (137) we deduce 

* = hW) + + -p- + —g- (141) 

from which, by virtue of Eq. (139), we obtain t (a) = 0. Now, with the aid of Eqs. 
(140) and (141), we may determine the components of the total energy-momentum 
tensor [Eq. (101)] attributed to a spacetime region due to the propagation of a null- 
string in it. As regards the corresponding (00) component, we have 

J°° = -Z I da J dr i 2 5(t - t') (142) 

In Eq. (142), we perform the substitution (r, a) — > [t(r, a), a] in a way such that 

dtda = idrda (143) 

thus obtaining 

J 00 = -Z J da Jdti 5(t - t') (144) 
Now, by virtue of Eq. (141) and the identity 

J dfi(k) f(k') S(k - k') = f(k) (145) 

Eq. (144) is finally written in the form 

T r 27T 1 / 

I 00 = T - / da — Jxl b 2 c 2 + yl c 2 a 2 + zl a 2 b 2 (146) 
2 Jo abc v 

In the same fashion, we obtain 



J 01 


= T 


T 
~2 


f 2w x (<t) 
Jo a? (t) 


da = 


jlO 


j02 


= T 


T 
2" 


f 27T yo{p) 

Jo b 2 {t) 


da = 


j20 


j03 


= T 


T 
2" 


f 2n z (a) 
Jo c 2 {t) 


da = 


j30 


/■27T 


c 






1 





(147) 



J 1 r 2n r 1 

I 12 = *o x oVo- h , o da = I 21 

2 JO ab . 1^2 tf. n 2 _|_ 1( 2 r 2^2 _|_ r 2 ^2^2 

T r 27r h 



2 Jo 



2 Jo 





b 2 c 2 


+ c 2 a 2 + zl 


a 2 6 2 






1 






b 2 C 2 


+ y 2 c 2 a 2 + z 2 


a 2 b 2 






1 





I 23 = Tl r //n :n f ' ^ l 31 (1 IS) 

" ' &c yx 2 6 2 c 2 + y 2 c 2 a 2 + *g a 2 6 2 
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and 



rii T f 27T 2 be 1 

2 « 3 yx 2 , 6 2 c 2 + yl c 2 a 2 + zl a 2 b 2 

122 = =f? / 1 ^ 

' ' &3 yx 2 , 6 2 c 2 + y 2 c 2 a 2 + 2 2 a 2 6 2 



2 Jo 



T r 2 ^ 9 ab 1 , . , 

J 33 = =F — / xg— flV 149 

2 ^ c 3 ^2 6 2 c 2 + y 2 c 2 a 2 + z 2 a 2 6 2 

Taking into account a circular (planar) string, namely 

x (a) = cosa , yo(a) = sina , z = const. (150) 

which propagates along the anisotropic z— direction of an axisymmetric Bianchi Type 
I model, i.e. 

a(t) = b(t) (151) 

we obtain t ^ t(a). Then, the only non-zero components of the corresponding total 
energy-momentum tensor are 

TFT 



J 00 = =f y° 2 + z o a2 

ac v 

jll _ ^ C _ j22 

2 « 3 . /c 2 + ^ 2 a 2 ~ 



7rT 2 a 

'lc 2 + zla? 



r : — : o , ' ,, , (152) 



and 



Z 03 = 2~ = / 3 ° (153) 

Due to the homogeneity of the spacetime under consideration, we may impose that 
Zq = [29]. In this case, the resulting total energy-momentum tensor is diagonal 

TFT 



r 00 = =p 



a(t) 



I 11 = =f — =/ 22 
+ 2a 3 (t) 

J 33 = (154) 

Now, the corresponding physical quantity to be inserted into the r.h.s. of the field 
equations, may be defined as: [T^ u ] = [1 1X11 per unit of proper — comoving volume], 
thus resulting to 



rp22 

(155) 



j>00 
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In this respect, the Einstein field equations which determine the evolution of a 
Bianchi Type I model with matter in the form of null-strings, are written in the form 

1 ^da^ 2 da 1 dc ^ .Gr. 1 (156) 

a 2 dt a dt c dt a' a 3 c 

1 d 2 a ^ 1 d 2 c ^ 1 da 1 dc ^ / G\ 1 (157) 

a dt 2 c dt 2 a dt c dt a' a 7 c 

2 I^ + I^I^ = (158) 
a dt z a dt c dt 

For the purpose of numerical analysis, in what follows we shall treat the dimensionless 
constant ^ (where a' is Regge slope) as a free parameter, to be varied at will. In fact, 
this parameter may be used to denote the regime at which gravitational phenomena 
dominate over the string effects (^ > 1) and vice versa (^ < 1). 

In principle, we may integrate the system of Eqs. (156) - (158) to obtain the exact 
form of the unknown scale functions, a(t) and c(t), thus determining the evolution of 
the cosmological model under consideration. Since this is not an easy task, we may get 
a good estimation of their dynamic behaviour through numerical integration. From 
Eqs. (156) - (158) only two are truly independent. The third one corresponds to 
an additional constraint to be satisfied by the solutions of this system. As such, we 
choose Eq. (156). The remaining independent field equations (157) and (159) may 
be recast in the form of a first order system, as follows 

= -Hl-\H lHi (159) 

and 

< ^ = aH l , d 4 = cH z (161) 

dt dt 3 K ' 

where Hi and H% are the anisotropic Hubble parameters. Accordingly, in what follows, 

we integrate numerically the system of Eqs. (159) - (161). 

To avoid possible implications of small variations on the dynamical parameters 
involved, the constraint (156) is checked to be satisfied with an accuracy of 1CT 10 
along numerical integration. Both Hubble parameters are measured in units of 10~ 4 
sec -1 , being normalized with respect to \[G. The initial conditions imposed on the 
dynamical variables of the problem (H®, H®, a , c ), are chosen so that: (a) a = c , 
i.e. at the origin, the two factor spaces (the isotropic xy-plane and the anisotropic 
^-direction) are separated, but of the same linear dimension, (b) Hf > 0, i.e. initially 
the isotropic space expands, in accordance to what we observe at the present epoch. 
As regards the anisotropic direction it may be either expanding or contracting, a thing 
that depends on the exact value of ^, i.e. on the regime at which we are working at 
(classical or quantum). 

The time coordinate is measured in dimensionless units, being normalized with 
respect to the Planck time, t — > t/t Pl (t Pl = \[G ~ lCT 43 sec). The limits of numerical 
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integration range from to — 10 to tf — 10 5 . The upper limit coincides with the origin 
of the GUTs epoch, tour — 10 5 tp;, corresponding to the end of the string era. On 
the other hand, the lower limit (to = lOtpz) is chosen as being safely away from the 
Planck epoch, since, in the absence of a self-consistent quantum theory of gravity, 
there is always a region of ambiguity around t — 0, of the order of the Planck time. 

The solutions to the system of differential equations (159) - (161) may be repre- 
sented as curves in a Hi — H% plane. Any point located on them, always satisfies the 
constraint condition (156), as well. Thus, these curves actually represent orbits of the 
dynamical system under study. Each curve, corresponding to a different set of initial 
conditions, is bounded by fixed points (or infinities) and represents a different type 
of evolution for the Universe. In what follows, we focus attention on the existence 
and the evolution of attracting points in the Hi — H 3 plane. The reason rests in the 
physical meaning of the attractor: No matter what the behaviour of a cosmological 
model at the origin might be, it will always end up to evolve as indicated by the 
location of the attracting point in the Hi — H 3 plane. 

In order to study a possible dissipation of the existing anisotropy, we furthemore 
define the so called anisotropy measure [30], as 

AH 3 Hi - H 3 . . 

This dimensionless quantity is a measure of the variance in the expansion rates at 
different directions and hence it represents the degree of anisotropy. Accordingly, we 
examine its dynamical behaviour during numerical integration. We expect that, at 
the end of the backreaction process, AH/H must settle down to a constant value 
that approaches zero, corresponding to a dynamical isotropization of the anisotropic 
background. 

Depending on the exact value of the free parameter ^, we consider three different 
cases: 



(i) Gravitational phenomena are comparable to string effects 

In this case, we have ^ = 1. The time evolution of the Hubble parameters for a 
particular set of initial conditions, (H®,H®) = (1,-1.5), is presented in Fig. 2a. 
The cosmic-time coordinate is measured in units of 10 3 tp/. We see that, although at 
the begining both Hubble parameters evolve anisotropically, after some time which 
is approximately of the order At £s 8 x 10 3 tpi, their evolution becomes identical. In 
the later stages (for t > 3 x 10 4 t P i), both evolution rates are reduced to a static value 
approximately equal to zero. In this respect, a complete isotropization of the model 
is achieved (i.e. Af = 0) quite rapidly (see Fig. 2b). Notice that the evolution rate 
of the anisotropic direction (H 3 ) never exceeds in value the corresponding parameter 
of the isotropic space (Hi). 

The above results indicate that, during the cosmological evolution there exists 
an attracting point in the Hi — H 3 plane. The explicit location of this point is pre- 
sented in Fig. 3. We observe that, for a wide range of initial conditions (including 
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both expansion and contraction of the anisotropic direction), every orbit of the dy- 
namical system under consideration ends up to a fixed point very close to (0,0), i.e. 
the Minkowski spacetime. This point is stable and according to the definition of 
the attractor, the cosmological model under study is not only isotropic, but also 
asymptotically (as t — > oo) flat. 

(ii) Gravity dominates over string effects 

In this case, we take ^ = 10. The time evolution of the Hubble parameters corre- 
sponding to the initial conditions of the previous case, is presented in Fig. 4a. The 
cosmic-time coordinate is also measured in units of 10 3 tp;. We see that at the end of 
the backreaction process, an initially large degree of anisotropy is subsequently min- 
imized, so that, asymptotically, the evolution rates of both factor spaces are reduced 
to the static value (0,0). However, in this case, the anisotropy dissipation rate is 
quite lower and the isotropization of the model is completed only after a period of 
At > 3 x 10 A tpi (see Fig. 4b). We furthermore observe that, although the dynamical 
behaviour of Hi is more or less the same as in the previous case, the anisotropic di- 
rection undergoes a period (~ 10 2 tpi) of accelerating (inflationary) expansion, during 
which the value of if 3 exceeds Hi. 

In this case also, we verify the existence of an attractor in the H\ — H 3 plane. 
The explicit location of this point is presented in Fig. 5. Notice that the dynami- 
cal behaviour of the various orbits is almost identical (there is no spreading of the 
orbits in the Hi — if 3 plane), although they correspond to quite different sets of 
initial conditions. Again, the attractor coincides to the stable point (0,0), i.e. the 
Minkowski spacetime. Therefore, the cosmological model under consideration is also 
asymptotically flat. 

(iii) String effects dominate over gravity 

In this case, we take ^ = 0.1. The corresponding results are more or less similar to the 
previous ones, with only one major difference: We do not have orbits corresponding 
to negative initial conditions of H%. For H® < 0, the constraint equation (156) 
is not satisfied any more. Accordingly, a dominant null-string structure does not 
permit contraction of the anisotropic dimension and therefore, the cosmological model 
under consideration does not admit a Kasner-type solution. Again, we verify the 
existence of an attractor in the H1—H3 plane, which coincides to the stable point (0, 0) 
corresponding to Minkowski spacetime, while anisotropy dissipation occurs rather 
rapidly (see Figs. 6 and 7). 

In concluding, we may state that the introduction of a null-string structure in 
a Bianchi Type I spacetime, may act in favour of both isotropization and flatness 
of the anisotropic background. Moreover, in the case where the string effects are 
either comparable or dominant with respect to the gravitational ones, the dynamical 
behaviour of H 3 indicates that the string actually decelerates the evolution of the 
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anisotropic direction (e.g. see Figs. 2a and 4a), thus serving as a factor for counter- 
inflation, as well (in connection see [26]). 

VII. Discussion and Conclusions 

In the first part of this paper the introduction of a generic ansatz has provided 
classes of exact, null-string solutions in a Bianchi Type I spacetime. This can be 
carried out irrespectively of the dynamical theory under which thecorresponding cos- 
mo logical solutions may originate and in fact, one has to use only the explicit func- 
tional form of the scale factors. This procedure is applied explicitly to the case of 
dilaton-Bianchi Type I cosmological models. 

We may realize the possibility of backreaction of the strings on the curved back- 
ground, by assuming that the dilaton and electromagnetic fields have the same cosmic- 
time dependence with the string. This specifies the values of the parameter /i, in such 
a way, that allows transition to a less or more anisotropic member of the same family 
of the cosmological solutions. The range of values of the cosmic-time parameter in 
which this process is valid, depends on the exact value of the string tension and the 
number of the string modes that are excited, in the first order approximation. If 
the three dimensional volume that encompasses the string is considered as a whole, 
its interaction with the spacetime geometry has a rate which is proportional to the 
parameter that characterizes the general class of solutions. 

Finally, if one assumes that in the primordial stages of the evolution of the Uni- 
verse, high spacetime curvature resulted in particle production [27], then one is led 
to examine if string-like properties of particles can also lead to anisotropy damping. 
In this respect, the strings, being considered as objects carrying stress and energy, 
seem to provide the means for anisotropy damping in epochs where the spacetime 
could have different expansion rates along different directions. This is examined in 
the last part of the paper, where the cosmological field equations regarding a null- 
string distribution propagating in a Bianchi type I spacetime, are integrated without 
any further assumptions. This integration indicates that, for a wide range of initial 
conditions, the null-string structure acts in favour of both the isotropization and the 
flatness of an, initially, curved anisotropic background. 

Acknowledgements: One of us (A.K) would like to thank the Greek State Schol- 
arships Foundation (I.K.Y), for financial support during this work. Financial support 
from the Greek Secretariat of Research and Technology under contract PENED-1768 
is gratefully acknowledged. 

References 

[1] M.Green, J.Schwarz, E.Witten, Superstring Theory. Cambridge University Press 
Cambridge 1987. 



23 



[2] H. J. de Vega Strings in Curved Spacetimes. Lectures delivered at the Erice 
School " String Quantum Gravity and Physics at the Planck Energy Scale" , June 
(1992), World Scientific. 

[3] H. J. de Vega and N. Sanchez String Theory in Cosmological Spacetimes, Lec- 
tures at the Erice School, September (1994). 

[4] A. L. Larsen and N. Sanchez Strings and Multi-Strings in Black Holes and Cos- 
mological Spacetimes in "New Developements in String Gravity and Physics at 
the Planck energy scale", World Scientific, (1995). 

[5] J. M. Salim, S. L. Sautu and R. Martins (1997), Class. Quantum. Grav. 15 1521. 

[6] D. Garfinkle, G. T. Horowitz and A. Strominger, Phys. Rev. D 43 3140 (1991) 

[7] G. W. Gibbons and K. Maeda, Nucl. Phys. B 298 741 (1988) 

[8] J. H. Home and G. T. Horowitz, Phys. Rev. D 46 1340 (1992) 

[9] R. Sorkin, Phys. Rev. Lett. 51 87 (1993) 

[10] D. J. Gross and M. J. Perry, Nucl. Phys. B 226 29 (1983) 

[11] A. G. Doroshkevich, Astrophys. J. 1 138 (1965) 

[12] K. S. Thorne, Astrophys. J. 148 51 (1967) 

[13] K. C. Jacobs, Astrophys. J. 155 379 (1969) 

[14] C. B. Collins, Commun. Math. Phys., 27 37 (1972) 

[15] P. I. Porfyriadis and D. B. Papadopoulos, Phys. Lett. B 417 27 (1998) 

[16] A. Schild, Phys. Rev. D 16 1722 (1977) 

[17] A. Karlhede, U. Lindstrom, Class. Quantum. Grav., 3 L33 (1986) 

[18] H. J. de Vega, A. Nicolaidis, Phys. Lett. B 295 214 (1992) 

[19] H. J. de Vega, I. Giannakis, A. Nicolaidis, Mod. Phys. Lett. A 10 2479 (1995) 

[20] C. O. Lousto, N. Sanchez, Phys. Rev. D 53 6842 (1996) 

[21] S. Kar, Phys. Rev. D 53 6842 (1996) 

[22] M. P. Dabrowski, A. L. Larsen, Phy. Rev. D 55 6409 (1997) 

[23] S. N. Roshchupkin, A. A. Zheltukhin, Class. Quantum. Grav. 12 2519 (1995) 

[24] C.W.Misner, K.S.Thorne, J.A. Wheeler, Gravitation. W.H. Freeman, (1973). 

24 



[25] F. Mandl, Statistical Physics, Wiley, New York, 1984. 

[26] A. A. Tseytlin, C. Vafa, Nucl. Phys. B 372 443-466 (1992). 

[27] N.D. Birrell and P.C.W. Davies, Quantum Fields in Curved Space. Cambridge, 
(1992). 

[28] J. V. Narlikar, Introduction to Cosmology, Jones and Bartlett Publishers, Boston 
(1983). 

[29] I. Giannakis, K. Kleidis, A. Kuiroukidis and D. B. Papadopoulos, Anisotropic 
Null-String Cosmologies. Accepted for publication in Mod. Phys. Lett. A (1998). 

[30] B. L. Hu and L. Parker, Phys. Rev. D 17, 933 (1978). 



25 



Figure Captions 



• Figure 1: Plot of the function f(p) — p(l —p)/(l + p) 2 

• Figure 2a: The time evolution of the anisotropic Hubble parameters Hi (squa- 
res) and i?3 (solid line), for a particular set of initial conditions [(Hf,H®) = 
(1, —1.5)], in the case where ^ = 1. The time coordinate is measured in units 
of 10 3 tpi. Notice that, the cosmological model becomes both isotropic (for 
t > 8 x 10 3 tpi) and flat (as t — >• oo). 

• Figure 2b: Plot of the anisotropy measure {^-) versus time (in units of 10 3 tpj), 
as regards the evolution of the Hubble parameters corresponding to the set 
of initial conditions (Hf,H®) = (1,-1.5), in the case where ^ = 1. The 
isotropization of the model is completed after 8 x 10 3 tp\. 

• Figure 3: The orbits of the dynamical system determined by the cosmological 
field equations, for five different sets of initial conditions, in the Hi — H% plane. 
All orbits approach asymptotically at the attracting point (0,0). 

• Figure 4a: The time-evolution of Hi (squares) and H 3 (solid line), for the 
same set of initial conditions as in Fig. 2a, in the case where ^ = 10. The 
time coordinate is measured in units of 10 3 tpi. Notice that in this case, the 
anisotropic direction undergoes a period of accelerating (inflationary) expansion 
before the model becomes both isotropic (for t > 3x 10 4 tp;) and flat (as t — > oo). 

• Figure 4b: Similar to Figure 2b, except that ^ = 10. Notice that the 
isotropization of the model is quite slower in this case. 

• Figure 5: Similar to Figure 3, except that ^ = 10. Again, all orbits end at 
the attracting point (0,0). 

• Figure 6a: The time evolution of Hi (squares) and H% (solid line) for a par- 
ticular set of initial conditions [(H®, H®) = (1,2)], in the case where ^ = 0.1. 
The time coordinate is measured in units of 10 3 tp\. Notice that the dominant 
null-string acts in favour of a counter-inflation of the anisotropic dimension, 
thus resulting to a rapid isotropization of the cosmological model. 

• Figure 6b: Plot of the anisotropy measure icjf) versus time (in units of 10 3 tp;), 
as regards the evolution of the Hubble parameters corresponding to the set of 
initial conditions H®) = (1, 2), in the case where ^ = 0.1. 

• Figure 7a: Similar to Figure 6a, except that (H®,H®) = (1,0). 

• Figure 7b: Similar to Figure 6b, except that (H?,H$) = (1,0). 
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